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Abstract 

We analyze the possibility of tomographic reconstruction of a sys- 
tem of three-level atoms in both non-degenerate and degenerate cases. 
In the non-degenerate case (when both transitions can be accessed inde- 
pendently) a complete reconstruction is possible. In the degenerate case 
(when both transitions are excited simultaneously) the complete recon- 
struction is achievable only for a single atom in the S configuration. For 
multiple S atoms, or even a single atom in the A configuration, only par- 
tial reconstruction is possible. Examples of one and two-atom cases are 
explicitly considered. 



1 Introduction. 

Quantum tomography is now an accepted technique to reconstruct the state of 
a quantum system [T] . Recent applications include the reconstruction of numer- 
ous physical systems, such as a radiation field [2], trapped ions and molecular 
vibrational states [3]- [5], spin [5] [7] and some other systems [5][H]. 

The main idea behind Quantum Tomography is to use population measure- 
ments of the "rotated" density matrix of the system. Explicitly, if p is the 
density matrix, a tomogram ui(tp, n) is defined by 

^,k)= (iP\U( K ) pU- 1 ^)^) (1) 

where 

P (k) = U(k) P U- 1 (k) (2) 

is the density matrix rotated by the unitary transformation U and k stands 
for all the parameters required to uniquely specify U. By varying the input 
parameters k, one obtains a complete collection of different observables (called 
a quorum), from which a characterization of the initial quantum state of the 
system can be obtained [TO] . 
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Every tomographic scheme is centered on the possibility of inverting Eq. (TT|) . 
This inversion and the corresponding reconstruction of p from u>(^>, k) can al- 
ways be mathematically achieved [TT] Q~2] when a Lie group (3 acts irrcducibly 
in the Hilbert space appropriate for the description of the quantum system 
under investigation. Experimental success depends on whether or not all the 
requirements element in (3 can be practically implemented. 

The objective of this paper is to study the possible tomographic reconstruc- 
tion of the density matrix describing a collection of three-level atoms. Several 
approaches have already been proposed to reconstruct the state of a general 
three-level system (qutrit) . In p2] Q3] [H] , the reconstruction of the quantum 
state of a three-level optical system is implemented for a frequency- and es- 
timating quantum states and measuring fourth-order field moments. The use 
of non-orthogonal measurements as a way to reconstruct the state of a system 
(provided those measurements span the Hilbert space) as well as a detailed 
example of reconstruction for one and two-qutrit systems, is considered in [7J . 

Even though one technically may require only a finite number of different 
experimental set-ups for complete tomographic reconstruction of atomic states, 
we will focus on the so-called redundant reconstruction, which implies a con- 
tinuous set of measurements "blanketing" all the parameter space. We justify 
this redundancy on the grounds that the reconstruction of a many-atom system 
would, in practice, require a large number of such discrete measurements. 

Our strategy is to investigate tomographic reconstruction of the atomic state 
by probing atoms through the application of a carefully selected sequence of 
dispersive and resonant electromagnetic pulses. 

Once an initial pulse of classical light has created a state of collective excita- 
tion in an ensemble of cold atoms, another pulse converts the atomic excitations 
into field excitations generating, for different atomic configurations, photons in 
a well-defined spatial and temporal mode |16j . The total number of such pho- 
tons is determined by photoelectric detection so that a probability of detecting 
a given number of photons can be directly related to a tomogram [TT] . 

In this work, our analysis will focus on two of three fundamentally different 
atomic configurations. They are the so-called A and S configurations, dis- 
tinguished by the presence of transition degeneracies from the generic non- 
degenerate configuration. 

In Sec[3] we discuss this non-degenerate case, where transition frequencies 
are essentially different for distinct atomic transitions and where each transition 
can be independently interrogated by a pulse of the appropriate frequency. In 
this case, the density matrix can be completely reconstructed. 

In Scc[4] and [5l the systems under consideration contain atomic transitions 
having the same frequency; it is not possible to interrogate every transition in- 
dividually. We will consider these cases at length and highlight the differences 
between the (global) symmetries pertinent to the description of these inequiv- 
alent degenerate atomic configuration. We will show, for these cases, that the 
density matrix cannot in general be completely reconstructed and that the par- 
tial information extracted from the measurements in the case of A and S atoms 
is essentially different. 
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2 u(3), su(3) and three— level atoms 

The Hamiltonian governing the evolution of a collection of A three-level atoms 
in a classical field has the form 

H = H + H 12 + H23 (3) 

where Ho is the free atomic Hamiltonian. and fly is the interaction term be- 
tween levels i and j. 

The terms in Eq.([3]) are most transparently analyzed by introducing a set 
{Sij]i,j = 1,2,3} of collective transition operators that satisfy the standard 
commutation relations of the u(3) algebra: 

[S^, Ski] = SjkSu — SuSkj ■ (4) 

Thus, the Hilbert space for our systems naturally decomposes into a sum of 
subspaces invariant under the action of the Lie group il(3). 

We further assume that the A atoms are indistinguishable and their states are 
fully symmetric under permutation of the particle indices. Hence, the possible 
states of our system belong to a single unitariy irreducible representation of 11(3) 
having dimension ^(A + 1){A + 2) and denoted by (A,0) in mathematics. The 
dimension of the Hilbert space is given by the number of ways of distributing 
A bosons in three modes. 

If we introduce the atomic basis 

{\nin 2 n 3 } ,ni,7i2,n 3 > , m + n 2 + n 3 = A} , (5) 

where rij denotes the population in the j-th atomic level, the matrix elements 
of S^ can be easily evaluated using the Schwinger realization: 

Sij i-> a \ a,j . (6) 

In the one-atom case, this yields 

4-li) = l^ (7) 

where the following identification 

|100>«-|1>, |010)«->|2), |001>H3>, (8) 

has been made. Throughout this work, we will assume the ordering E\ < E2 < 
E3 of individual atomic levels. 

In terms of 5V/'s, the Hamiltonian of Eq.Q takes the form 

3 

(9) 

(10) 
(11) 





— 2^ ^ ^ u ' 
1=1 




H12 


= 5i(e lWlt £ 12 - 


h e- luJlt S; 


H23 


- .9 2 (e ?w £ 2 3- 


h e- luJ2t S : 
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where cJi and 0J2 are frequencies of the external fields and g\ and g 2 are coupling 
constants, chosen to be real for simplicity. 
The operator 

3 

= J2 S u (12) 

i=i 

commutes with all other operators in the u(3) algebra. This operator is pro- 
portional to the unit operator when acting on occupational states of the form 
\nin 2 n 3 ), 

N |nin 2 n 3 ) = A|riin 2 n 3 ) . (13) 

Removing N reduces u(3) to su(3). Thus, the possible evolutions generated 
by the Hamiltonian H are, up to an unimportant global phase, finite Sil(3) 
transformations. 



3 Non-Degenerate Case 

The non-degeneracy condition is understood to imply that the atomic transition 
frequencies (E3 — E 2 ) and (E% — E%) arc sufficiently distinct to satisfy 

(E 3 - E 2 ) - (E 2 - E ± ) > 9l , (E 3 - E 2 ) - {E 2 - E x ) > g 2 . (14) 

In this manner, each transition can be interrogated separately by an external 
field. A typical non degenerate system is illustrated in FigJTJ 




Figure 1: A typical level scheme for the non-degenerate case. 

For tomographic purposes, the level pairs (1 <-» 2), (2 <-» 3) and their re- 
spective transitions, if taken in respective isolation, could be considered as in- 
dependent two-level subsystems. However the full system must be treated as a 
three-level system under 611(3) evolution. 

In the rotating frame, the Hamiltonian of Eq.([3]) takes the form 

Hint = A23S33 - A12S11 + 312 (§12 + S21) + .g 23 (S23 + 532) , (15) 
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where 



A12 = E 2 — Ei — u>i 



A23 = E 3 



E 2 



■ U> 2 



(16) 



Because the frequencies uj\ and oj 2 of the external field are adjustable pa- 
rameters, two types of field pulses can be applied to our system. The first type 
is characterized by Ay « and thus resonant; it stimulates the corresponding 
atomic transitions. The second type is characterized by Ay 3> gij and is thus 
far-off resonant (dispersive); this kind of pulse leads to some phase shifts. 

The corresponding resonant and dispersive evolution operators have the form 



uR (^12) 

33 (023 ) 



ifiij (Sij + Sji 



UP 



cxp 
exp 

exp 1023533 , 



(17) 
(18) 
(19) 



where Pij = g%jtij and 0y = (Ay + g^/Ay)^. Here, Uj,Tij are time intervals, 
not necessarily equal. It should be noted, that for short interaction times Ty- 
satisfying g 2 nj/Aij « 1 and g 2 << Ay, the second term in the expres- 
sion for 



exp(— iAij iij 



ij can be obviously neglected. However, for long interaction times, 
Tij) becomes strongly oscillating and the measurements should be 
carried out in stroboscopic times, Ty = 27m/ Ay. 

For a complete reconstruction of the density matrix in the absence of degen- 
eracy, it suffices to measure the probability of detecting zero photons (i.e. zero 
fluorescence condition) in the irradiated field. This corresponds to the detec- 
tion the atoms in the ground state, |A00). It is worth noting here that, in the 
non-redundant scheme when only a finite number of different pulses have to be 
applied, the measurement of non-zero photons are required [17] . 

It is possible, by combining operations in Eqs. (TTT|) - (rr9")) . to obtain a element 
of the group il(3) sufficiently general for our purpose. This is best seen by 
first observing that an element of Sil(3) is parameterized by eight real num- 
bers and can be conveniently factorized [TS] into a product of 611(2) subgroup 
transformations : 

j7(ai,/3i,7i, a 2 ,P 2 , a 3 , p 3 , 73) 

= #23(0*1, 71) • i?i2(a 2 ,^2,a 2 ) • i?23(«3, $3,73) 7 (20) 

The action of such a group element on basis states of an irreducible representa- 
tion is given in [TB] . 
The notation 



n=(n 1 n 2 n 3 ), I = J 23 = \{n 2 + n 3 ) , 



(21) 



will be a useful shorthand throughout this paper. In particular, we note that 

(u\R 23 (a,p,j)\n) = {viv 2 v 3 1 R^{oi, P, 7) \nin 2 n 3 ) = ^m'.mJ^) > 
(v\R l2 (a,p,j)\n) =(^2U 3 \R l2 (a,p,j)\nin 2 n 3 ) =V I n \l„ l jn), (22) 
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where V J MM , an 611(2) Wi gner 23-function and 

I12 = §("i + n 2 ) , M„ = i(n 2 -n 3 ), m„ = |(ni - n 2 ) . (23) 

More generally, the definition of states for the irrep (A, /i) and 611(3) elements 
between those states can be found in [18]. Specialized results are collected inlAl 
In particular, some formulae in this section implicitly depend on matrix elements 
of states in irreps of the type (^4, 0), (0, A) and (cr, a). The notation of 611(3) T> 
functions conforms to that of [18] and uses the pair (A, fi) with the labels n, I 
to unambiguously distinguish the 611(3) T> functions. For occupational states, 
n, I are given by Eq. (j2"T]) . 

Because the state |A00) is, up to a global phase, unchanged by the action 
of dispersive pulses and operations of the form i? 23 , it is enough to consider a 
sequence of pulses of the form: 



E%23,#!3,¥>i2,/?i2) = U^i-^UH^UiiMU^n) . (24) 



In the single atom case, the evolution operator Eq. (|2^|) is a 3 x 3 matrix 
explicitly given by 

U{f23,p23, ¥>12, A2) 



U(^f23, /?23, ¥?i2) P12) can be more easily analyzed in the factorized form 



/ 



e^ 12 cos (/3i 2 ) 
-i cos (flzs) sin (/3 12 ) 
e" 1 ^ sin (/3 u )sm (/3 23 ) 



-ie^ 12 sin(/3i 2 ) 
cos (J3u) cos (/3 23 ) -« sin (/? 23 ) 

ie"^ cos (/?i 2 ) sin (/3 23 ) e" 4 *^ cos (/? 23 ) 



V 



(26) 



where 



U 




and 



X = :j(2</?i 2 + ^23) ■ 



(28) 



Clearly, the matrix C/(</? 23 , /3 2 3, v?i 2 , /3i 2 ) of Eq. (|2"6")) is an element of 11(3) whereas 
U of Eq. (|27p is an 611(3) transformation. Comparing with the paramctrization 



of [TS], we have the correspondences 

OLl -> -1^23 — §7T , /3i -> 2/3 23 , 71 + |^12 + §V>23 

OL2 —> -§¥>12 - |v?23 , #2 -> 2/3i2 , 

Q!3 -3^12 - 3V23 - 7T , 03 = 0, 73 = 0. 

(29) 

Note that, although U is not the most general 6iX(3), it can be multiplied on 
the right by a transformation of the form 

i?12(a2,0,a 2 )i?23(a3,/?3,73) = 

e~ ia2 \ 

e ii(2a 2 -a 3 -73) cos (^) _ e ^(2a 2 -a 3 +73) sm (|i) (30) 

e i(&3-^) sin(&) e5( S3+ ^)cos(%) / 

without affecting the dynamics of the 1 100) state. Thus, U is equivalent to a 
general transformation when acting on 1 100) . 

Expanding the density matrix in the occupational basis: 

P = ^2 \ n )( v \Pn,», (31) 

1 1 V 

and introducing the shorthand r — {^iz^iz^n^vi) > we rapidly obtain 

= E(- 1 )" & o ) )o >n zW^( ( o°k A i)o ) ,.Z'WPn,,, (32) 

where 

P £'W 2 ( T ) = ((A,/i)n 1 J 1 |Cr(T)|(A,/i)n a J 2 > (33) 

is an 6il(3) Wigner-function for the irrep (A, [i). Further notational details and 
properties of these functions (in particular Eq. ([7D|) ) can be found inlAl 

Products of <5il(2) P-functions can be decomposed into sums of P-functions 
multiplied by products of 6H(2) Clebsch-Gordan coefficients. The same holds 
for products of 6il(3) P-functions provided that we use Sil(3) Clebsch-Gordan 
technology. Thus, given that 6il(3)-coupling (A, 0) <g> (0, A) decomposes in the 
direct sum [T^j 

(A,0)®(0,A) = (A, A) ®(A _ 1,^-1) ©...©(0,0), 

A 

= ©(A, A), (34) 

A=0 

we have 



V 



(A00)0,n/\' I ^(0AA)0, v*I'\' > 
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v-WU) , , / (AO) (0,A) ^a,aj v 

2^ (aaa)o,jvWjV^ ^(^00)0 (oAA)O (AAA)O / 



(A, A) \ 



(AO) (0,A) 



(A, A) 



nl v*I' 
where 

iV (A) = (m + i/J - (A - A), n 2 + i£ - (A - A), n 3 + ^3 - (A - A)) . 
and where 



(AO) (o,A 



(A, A) 
NWl 3 



(35) 



(36) 



(37) 



is the 611(3) Clcbsch-Gordan coefficient for the coupling of \(A,0)niIi) and 
|(0, A) 712/2) to I (A, A) N^'Ia). The appearance of extra factors of (A — A) in 
the construction of ' is discussed in Eq. ([86|) of IBl 
Inserting Eq.flHU in Eq.(|32|) yields 



w(r) 



2 (-ir^2? 



(A.A) 

(AAA)O, JV(A) J 



(r) 



711/ A J 



/ (AO) (o,A 
\(A00)0 (oaao 



(A, A) \ /(AO) (o,A 
(AAA)0 / \ 7i/ v*! 1 



(A, A) 



.(38) 



After some straightforward manipulations detailed in [Cj we obtain the final 
expression 



(/x + i) 3 / (AO) (o,A 

1024tt 5 \(A00)0 (0AA)0 



// .7 

(AO) (o,A 

77/ V*V 



(m>m) 
(w) 



^^ 0| „ w >) W (r).(39) 



As there is no restriction on n or 1/*, Eq. (|39p shows that, in the non- 
degenerate case, the density matrix can be completely reconstructed. 



4 Degenerate A-type atomic systems 

Let us turn our attention to the case of a degenerate A-type system. A typical 
A atom is schematically illustrated in FigJ^l 

In the single-atom case, the allowed transitions are |1) «-> |3) , |2) «-» |3) . 
The degeneracy condition is 

E 3 -E 1 = E 3 - E 2 . (40) 

In the multi-atom case, the only atomic configuration that can be unam- 
biguously identified by photon counting is when every atom is excited, i.e. when 
the system of A atoms is in the state |00A- 
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Figure 2: Schematic representation of a single atom in a A configuration. 
4.1 The evolution 

In the rotating frame, the interaction Hamiltonian has the form 

H A = AS 33 + g (S 13 + S 23 ) + g (S31 + S 32 ) , (41) 

where A = E 3 — E\ — oj and g\ = 92 = g for simplicity. 

In the single atom case, H A can be represented as the following 3x3 matrix: 



H A = g . (42) 




A simple basis transformation 

|1> - |1) = -^(|1)-|2)), |2)^|2) = -L(|1) + |2)) > 
1 3) ~ |3), (43) 
given by the constant matrix 

73 7f 

r « = I -7f 7f I - 


transforms Eq. l|32"j) to the block diagonal form 

/ \ 
H A ^H T = Tf/ ff A Ti 2 = V2g \ . (45) 

\0 V2.9 A / 

The effect T12 on basis states is illustrated in Fig|5J T\i produces a dark state 
|1) completely decoupled from the remaining doublet. In view of this we can 
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I3> 



ll> I2> 

Figure 3: The basis states resulting from the transformation T\%. 

expect, on general grounds, that a complete reconstruction will not be possible 
as our Hamiltonian H\ cannot possibly probe the dark state 1 1) . 

Using the basis {|1), |2), |3)}, the resonant pulses, with A = 0, are of the 
form 

/ 1 \ 

U£(V2gt)=\ coa(V2gt) ism(V2gt) . (46) 
\ isai(y/2gt) cos(V2gt) J 

In the same basis, the dispersive pulses, with A » j, are described in the 
stroboscopic approximation by the effective evolution operator 

1 

(2. 9 2 i/A) = | e 2 ^ 2 */ A (47) 

e ~^9 2 t/A 

In the two-dimensional subspace spanned by |2) and |3), the operators 
and C/jP correspond to Gil(2) rotations about the x and z axes, respectively: 

U*(a) ^ B% 3 (a) , UZ(P)^RW), (48) 

in an obvious notation. A sufficiently general sequence of pulses can thus be 
written 

U A (Cl) = ^(a). ^3(^-^3(7), 

= i?^ 3 (a + |) • i?f 3 (/3) • i?| 3 ( 7 - f ) = R 23 (fl) . (49) 

For the one-atom case, the 3x3 matrix representation for this evolution has 
the form 

- - ( 1 ° °\ 

C/ A (fi) = * * (50) 
\ * * J 

where * indicates a non-zero entry. The block diagonal form of U is explicit. It 
shows that, for a system containing one or more than one atom, there will always 
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be at least one subspace which cannot be reached in the course of evolution of 
|0(L4); such decoupled subspaces are an obstruction to complete reconstruction. 

To illustrate this, we expand the density matrix for an A-atom system in 
the basis {\n1h2n3)} of occupation of the states |1) , |2) and |3): 

P= Inifo^ih^hlp&MnaliihWa) • ( 51 ) 

hih 2 n 3 C'i v 2 v 3 

Using the shorthand n for the triplet hi-fis, we can write the tomogram as 
w(fi) = ^((XL4|fl 23 (f2) |nin2^3>(^i^2^ 3 |i?23(^)|00A) p^h^)^^) • (52) 

As the i?23(0) rotation does not affect the first atomic index and acts irre- 
ducibly as an 611(2) rotation in each the subspace spanned by states having a 
common fixed n\ , we write 

\nin 2 n 3 ) \h\\Im) , | v\Vi^z) -* | v\;Ip) , (53) 

where 

1= |(n 2 + na) = \{u 2 + v 3 ), m = \{n 2 - n 3 ) , fx = \(yi - v 3 ) (54) 

Furthermore, the tomograms of Eq. ([52|) must have h\ = V\ = 0, so Eq. ([54ll 
leads to 

m (i 

= ~Y] (-l)~ J ~ M Cf 1,1 -I C I m, I -u^OM^) P(0;Im),(Q;IH,) ' ( 55 ) 
m /J.L 

where P(„ 1MW ) -* P( ni ;im),( vi ;in) has been used to conform to the 

notation of Eq. (|53"l) , where 2?^, m / is the usual 6il(2) V function and where 
C L^ mi ,L 2m o is an eiI ( 2 ) Clebsch-Gordan coefficient. 

Multiplying both sides of ([55]) by 2?q ^-(fi), integrating over (511(2) and using 
orthogonality of the Clebsch-Gordan coefficients rapidly gives the elements of 
the density matrix that can be reconstructed from the tomographic process as 

~ _ / 1 \i+u \ "* 21/ - 1 „ LM (r<L0 \~ 1 

P(0;/m),(0;/A») — / y g_ 9 °Jra, J -/t 1°/ J, J -I) 

x y (56) 

4.2 Reconstruction for state of one and two A-type atoms. 

In a system containing a single atom, the Hilbcrt space is spanned, in the 
notation of Eq. ((53|) . by states of the form \n\ \ Im), with 

|1) = |100), |2) = |0;i|), |3) = |0;i,-i). (57) 
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Under the evolution Ua(Q) = #23(0), the initial state |0; \,—\) cannot reach 
the dark state so it is only possible to reconstruct element of p of the form 

P(0;^m),{a;-m')i w ^ tn m ' m ' = ^ T ^he ^ as * diagonal element, P(i ; oo),(i ; oo) can 
be inferred from the normalization. None of the remaining four matrix elements 
can be determined by our scheme. 

In the two-atom case, an even smaller proportion of matrix elements can 
be recovered. Using again the notation of Eq. (|53|) . states of the irrcp (2,0) 
are conveniently given, in the occupational, tensor product and (511(2) basis 
\ni;£m), in table [TJ 



Table 1: Basis states for two atoms in the A configuration. 



\nin 2 n 3 ) 


|l) ni |2)" 2 |3)" 3 


\hi] I, m) 


[002) 


|3)|3) 


[0; 1,-1) 


Oil) 


^(|2)|3) + |3)|2)) 


1 0; 1,0) 


|020) 


|2)|2) 


|0;1,1) 


101) 


^(|1)|3) + |3)|1)) 


I1-- --) 


|110) 


^(|1)|2) + |2)|1)) 


1 l- i I) 

1 ' 2 ' 2 / 


[200) 




1 2; 0, 0) 



The initial state |002) will not evolve out of the I — 1 subspacc, so only 
matrix elements of the form /)(o ; im)(0;im') can be reconstructed using our scheme. 
These represent only nine of the possible 36 elements of the density matrix. 

We conclude this section by noting that the situation obviously worsens (in 
the sense that a smaller and smaller proportions of the matrix elements can be 
recovered) as the number of atoms increases. 

5 Degenerate S-type atomic systems 

Finally, we consider the case of the S system. It is illustrated, for a single atom, 
in FigfJ] For this configuration, the condition E2 — E\ = E3 — E2 holds. 

5.1 The evolution 

In the rotating frame, the Hamiltonian governing the evolution of a collection of 
A atoms in the 5 configuration in an external field has the form (g\ = 92 = g) 

H3 = A (533 - S n ) + 9 (S 12 + S 32 + S21 + S 23 ) , (58) 

where A = (E 3 - E x )/2 - lu. 

Important insight into the nature of this Hamiltonian can be gained by 
noting that the operators Sn — S33 and S12 + 5*21 + "S23 + S32 are, in fact, 
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I3> 



2> 



ll> 

Figure 4: The S configuration for a single atom. 

proportional to two of the three generators of the so (3) subalgebra of su(3) : 

Sn - S33 !->■ L z , §12 + S 2 i + S23 + S32 i-> V2L X . (59) 

Thus, the possible evolutions are elements of the SD(3) subgroup of 611(3). 
Clearly, a convenient sequence of pulse is given by 

U S (Q) = R z (a)-R x (f3)-R z ( 1 ), 

= R z (a + ir/2)-R v (!3)-R z (~f~ir/2). (60) 

Here, the resonant pulses are of the form 

R X (P) = cxp (5i2 + S21 + S 23 + S 32 )/V2) (61) 

while the dispersive pulses are generated by S%i — 6*33. 

It is important to note that, in the one-atom case, the exponentiation of H-e 
in Eq. (|58[) produces an evolution that acts irrcducibly on the non-degenerate 
states of the Hilbert space: in contrast with Eq. (|4"9")l of the A case, the "ro- 
tations" R x and R z of the 2 states are not restricted to a two-dimensional 
subspace of the whole Hilbert space. 

To analyze the many-atom case, we start by observing that the state of the 
system for which every the atom is completely excited, |0CL4), is an eigenstate of 

L z with eigenvalue — A and is annihilated by Z_. Here, L_ = ^S*2i + S32) /V2 

is constructed in the usual way: L_ = L x — iL y . Thus, |00^4) is the unique 
angular momentum state \L, —L), with L = A: 

\00A) -> \L = A,M = A) . (62) 

As this state contains the largest possible number of excitations, it can be 
uniquely identified through photon counting so that the corresponding tomo- 
gram is determined from the probability of detecting 2A photons in the irradi- 
ated field. 
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The general correspondence between the occupational basis states In.1n.2n3} 
is found in [20] and given by 



\LM) = 



2 L + M (\{A + L))\{L + M)\{L - M)\(2L + 1) 
' (\{A- L))\{A + L + 1)\ 

U2 9 1 ty^- L y (a\r(4) L+M ~ 2p (4) p - M , n w fi o, 

S al) ~ 2a ^) ^ 2P p!( D -M)!(L + M- P )! |0>(63) 



Table 2: Angular momentum basis states as linear combinations of occupational 
number states for one and two atoms in the S configuration. 



A 


L 


M 


\LM) 


1 


1 


1 


al|0) 









410) 






-1 


4|o) 


2 


2 


2 


^(4) 2 |o) 






1 


44 10) 









1 

v/3 


((4) 2 + 44) 


10) 






-1 


44 10) 






-2 


^(4) 2 lo) 


2 










'(4) 2 - 244 


)|0) 



It is clear that, given an angular momentum state in the irrep (A, 0) of su(3), 
we can unambiguously write it as a linear combination of occupational states, 
and vice versa. Thus, we may expand 

P = \LiM 1 )(L 2 M 2 \p LlMl ,L 2 M 2 , (64) 

L 1 M 1 L 2 M 2 

where L\, L 2 run from A, A — 2, 1 or depending if A is even or odd. 

As the evolution is necessarily an element of 6:0(3), the tomogram takes 
the general form 

w(fi) = ]T (A,-A\R(n)\L 1 M 1 }(L 2 M 2 \rt{n)\A,-A}p LlMuL2 M 2 , 

L 1 M 1 L 2 M 2 

= ^ (^) I3 -1,M 2 (^) PAM U AM 2 , 

Mi M 2 

= ^2 {-iy A ~ Ah T>^ M (n)C^°_ Ai AA C AM U AM 2 PAM lt AM 2 , (65) 

with Ci 

l 2 m 2 a regular angular momentum Clebsch-Gordan coefficient. 
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In a manner similar to the previous cases, multiplication by Dq A ^(f2), inte- 
gration over SD(3) and orthogonality of Clebsch-Gordan coefficients yields 

i i \A+M 2 „ 

.7 J 

The result clearly shows that only those linear combinations of occupational 
state that transform by angular momentum L = A can be reconstructed. 

5.2 Examples: one and two atom cases 

For a single atom, we see from table[2j that the tomogram is constructed from an 
L = 1 state. There is no other angular momentum multiplet and so the evolu- 
tion, an clement of 6D(3), will yield sufficiently many tomograms to guarantee 
complete reconstruction. 

The matter is different for the two-atom case. In this case, the tomogram is 
constructed using an L = 2 state but the Hilbert space also contains an L = 
subspace, which cannot be reached from L = 2 with our evolution. Thus, if we 
are limited to measuring a total of 2A photons, it will only be possible to recover 
p2M,2M' and impossible to reconstruct poo,2M, Poo.oo, P2M,oo- This is because, 
in our scheme, it is not possible to extract photons from the L = state, and 
absence of photon does not pin down a particular state. 

It may be possible to measure fewer than 2A photons, but this does not 
lead to more information. There is only one state with M = A and one state 
with M = A - 1 (or M = -A and M = -(A - 1)). It is possible to use the 
L = A, M = A— 1 state for the tomogram and measure 2A— 1 or 2 A— 2 photons, 
but we will recover nothing more than if we had started with L = A, M = A. 

There are two states with M = A — 2; they belong different angular mo- 
mentum multiplet. Thus, if we measure, say, a total of 2{A — 2) photons, it is 
not possible to know unambiguously if this is the result of a complete cascade 
within the L = A — 1 multiplet or a partial cascade within the L = A multiplet. 
This kind of limitation becomes obviously more severe as the number of angular 
momentum multiplet containing a given M value increases. 

6 Conclusions 

We have proposed a physical realization applicable to the reconstruction of the 
quantum state of three-level atomic systems. The information about atomic 
states is extracted by measuring the total number of excitations after successive 
applications of electromagnetic field pulses. 

We have shown that, in the non-degenerate case, the complete reconstruction 
of atomic states is possible. Although the number of independent parameters 
required for a complete reconstruction is less than needed for the complete 
parametrization of a generic clement of 611(3) group, a complete reconstruction 
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is possible because, in addition to the usual evolution of the system, another 
tool is available in the reconstruction scheme: the projective measurement. 

When degeneracies are present, the possibilities of reconstruction are lim- 
ited. The origin of these limitations is essentially different for atoms in A and 3 
configuration. In both cases, the evolution operator operators belong to a sub- 
group of the whole SiX(3) group, and our work illuminates the subtle distinction 
between the global properties of (5,0(3) and @il(2) as subgroups of SH(3). 

In the S case, the reconstruction is rooted in an 60(3) symmetry of the phys- 
ically available evolution operator; this symmetry provides information about 
a single subspace. In the one-atom case, the Hilbert space contains precisely a 
single 60(3) subspace, so the density matrix can be completely reconstructed. 
In the multiple-atom case, only reconstruction in one pre-determined subspace 
is possible. In this case, our protocol would be to apply the sequence of pulses of 
Eq. (|60|) with a subsequent measurement of the number of 2 A of photons in the 
irradiated field, giving us the tomogram appearing in the reconstruction formula 
of Eq.(jnH). For completeness, we note here that we did not consider the effective 
two-photon-like transition in the 5 system due to extremely narrow width of 
such transitions (~ g 2 / A), which leads to serious experimental difficulties in its 
detection. 

In the case of A configuration, the evolution operator generates an ©it(2) 
transformation and, even in the one-atom case, there is always more than a 
single Sil(2) multiplet: a complete reconstruction is impossible because there 
always exists invariant Sil(2) "dark" subspaces, which cannot be uniquely iden- 
tified by measuring irradiated photons. We stress that the decomposition of the 
Hilbert space into invariant subspaces occurs as a result of the inability to access 
independent transitions separately; this to be contrasted with the approach of 
Ref.[S], wherein 6il(2) decomposability arises from considerations of perfectly 
general polarization states. Note also that although the effective transitions 
between degenerate levels in the A case are not sensitive to the atom-field de- 
tunings, they still require long interaction times. 

The tomographic protocol for A differs from the 2. After application of the 
sequence of pulses Eq. (f2~4")) to a A-type atom, we have to measure the probabil- 
ity of detecting zero irradiated photons, which leads to the tomogram used in 
Eq.©. 

Finally, we observe that the tomographic reconstruction process for a col- 
lection of non-degenerate three-level atoms is a simple generalization of the 
familiar process used for two-level quantum systems. In both instances, one 
uses the whole dynamic symmetry group to carry out the inversion process. 
In contrast to this, we are restricted to a specific subgroup in the degenerate 
cases, which essentially reduces our tools and actually limits the possibility of 
the complete tomographic reconstruction. 

We would like to thank Dr. O. Aguilar for his participation in the early 
stages of this project. The work of A. B. Klimov is partially supported by 
grants CONACyT 45704. The work of H. de Guise is supported by NSERC of 
Canada. 
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A (511(3) basis states and D-functions 



In this section, we review some notation useful mostly in section [3] Further 
details can be found in [TH] . 

The correspondence between the occupational basis and states of the (A, 0) 

is 

|nin 2 n 3 ) >-> |(4,0)n/>. (67) 

In Eq.(|67|) and throughout this paper, n is a shorthand for [jixn^n-i). Here, the 
(A, 0) labels indicate that \n1n2n3} can be reached, using the S!y operators of 
Eq.©, from the state |j400) . This state is killed by the so-called su(3) raising 
operators S±2,Sx3 and 523. The eigenvalues of the su(3) diagonal operators 

f)i = Six — S22 , t)2 = S22 — S33 , (68) 

acting on \A00) are, respectively, (A, 0). 

The angular momentum label / is necessary to deal with the general case 
considered in Ref. [18j , where states of more general families of the type (p, q) 
are constructed. A state \(p,q)nl) can be reached from the state \(p,q)(p + 
q,q,0)jp), i-e. with m =p + q,n 2 = q,n 3 = and / = \p. \(p,q)(p + q 1 q,Q)\p) 
is killed by the su(3) raising operators, and the eigenvalues of (hi, h 2 ) arc (p, q). 
When p and q are both non-zero, it is possible to have distinct states in the 
same (p, q) family that have identical rt, so the angular momentum label I is 
required to distinguish these distinct states. 

Some calculations require the evaluation of the matrix elements 

(nxn 2 n 3 \U(a)\nxn 2 n 3 }* - {(A, 0)nl\ U (a) | (A, Q)vl'}* , 

= (69) 

This matrix element is related to the matrix element between basis states of the 
irrep (0, A), which is conjugate to (A, 0), by 

= (-1)" 2+ ^ V { XlAr) . (70) 

Here, 

V^ v , r {r) = ((0,A)n*I\U(a)\(0,A)v*I>) . (71) 

The dimension of (0, A) is the same as the dimension of (A, 0), but the con- 
struction of Ref. [18] for basis state of |(0, A)n*I) requires twice as many quanta 
as the basis states of \(A,0)nI). The relation between n and n* is 

n = (711,712, n 3 ) t-* n* = (A - rix, A - n 2 , A - n 3 ) . (72) 



Using this and the results from [18], one can verify Eq. ([70|) . 

Using Eq. ([70]) . one can also verify that the T> functions are orthogonal, in 
the sense that 

f dn©^(n)2?^ L ,(n) = ,y AV v*nn'^^^, (73) 
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where 

dfl = sin pi cos \ (sin 3 sin /3 3 da.\ d(3\ d'ji da 2 dfii da 3 <% CZ73 (74) 

. The nor- 



is the invariant measure, which can be found in the usual ways p 
malization follows from the dimensionality formula 

dim(A,//) = ±(A + l)0 + l)(A + /i + 2) 

for the irrep (A, [i) and the use the parameter range 

< ai < 4tt, < /3i < 7r, < 71 

< a 2 < 2tt, < /3 2 < 7r, 

< a 3 < 4tt, < /3 3 < 7T, < 73 



< 



< 



(75) 



47T . 
47T . 



(76) 



B Reduced (3il(3) Clebsch-Gordan coefficients 

B.l Basis states 

The construction of states in the irrep (p, q) of su(3) is detailed in [18] . We can 
summarize this procedure by stating that one requires, at a minimum, a total 
of p + 2q bosons. These bosons must be of at least two types when q ^ 0. Thus, 
if a\, creates a boson of type j in mode i, we define, quite generally, 



,t 



,t 



Ski = a kl a H + a k2 ae2 ■ 

If |0) denotes state with no boson excitation, the state 



11 
'21 



12 
'22 



(a t 11 ) p |0)~|(p, g )(p + g ,g,0);ig) 



(77) 



(78) 



containing p + q bosons in mode 1, q boson in mode 2 and none in mode 3, 
belongs to the (p, q) irrep. It is, in fact, killed by every Cki with i > k and is 
thus the highest weight state of (p, q) . Here, 



'11 



'12 



u ll u 22 



a 1 a} 

u 12 u 21 



(79) 



is the determinant of the matrix. Other states in (p, q) are obtained by laddering 
down from \ {p, q)(p + q, q, 0); \q). 

This is not the only possibility. One can verify that 

\(p, q)(p + q + t,q + t,t);hq) 



"21 

4i 

°21 

4i 



'22 

':)2 



"12 

'22 

,t 
'32 



' 2:5 
'•33 



'13 

,t 
'23 



'21 



'22 



(aliFlO), 



l(Pi q){p + q,q,oyAq) 



(80) 
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visibly contains p + 2q + 3i bosons but is equivalent to the state of Eg. ([75)1 
because the determinant 



,t 



'13 



"21 "22 "23 

Cto-i (loo (In 



(81) 



-"31 "32 "33 

is an 611(3) scalar. 

More generally, if the usual ket \(p, q)nl) denotes a state in (p, q) containing 
p + 2q bosons, then the (round) ket 



\{p, q)nl) 



4i 

4i 



'22 

,t 

'32 



'13 

,t 
'23 



'33 



\(p, q)nl) 



(82) 



differs from |(p, q)nl) by at most a normalization but contains p+2g+3i bosons. 
B.2 (5il(3) Clebsch-Gordan coefficients 

The Sil(3)-coupling (A, 0) (8) (0, ^4) can be decomposed in the direct sum [19] 
(.4, 0)8(0, A) = (A,A)®{A-1,A-1)®...(B(0,0), 

A 

= ©(A, A). (83) 

A=0 

The irrcp (a, a) occurs at most once in the decomposition. 

To compute (511(3) Clebsch-Gordan coefficients for states in the series of 
Eq. ([53"[) . we must couple states of the form 



|(A 0)nIi)\(p,A)v*I 2 ), 



(84) 



which contain a total of 3v4 bosons of three types. States in the irrep (a, a) of 
the series of Eq. |83[) are of the form 



(a,a)NI 3 



'ii 



'31 



,t 
'12 



'32 



''13 
'23 
l 33 



A-c 



\{a,a)N^h) 



(85) 



where | (cr, a)N^I^) is the state with 3a bosons described in [18]. 

Note that, because the state |(cr, a)N^I^) does not contain 3A bosons, we 
do not have ni + v* = iV,- etc but rather 



N 



(<0 



(A- a) 



(86) 



Thus we have 



(a,a)N^I 3 ) = Yl l(A0)nJi)|(0,i4)t/*/ a > 



(AO) (0,A) 
nil v*I 2 



(a, a) 



(87) 
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Table 3: 1-atom case. Reduced 611(3) Clesbsh-Gordan coefficients for (1,0) ® 
(0,1) -(0,0). 



Ni 


h 


n\ 


h 




h 


/(l.O) (0,1) 


(0,0) \ 











i 

2 


1 


1 

2 


-v 


I 






1 











+v 


1 



Table 4: 1-atom case. Reduced 611(3) Clesbsh-Gordan coefficients for (1, 0) £g> 
(0,1) -(1,1). 



Ni 


^3 


n\ 


h 




h 


/(l.O) (0,1) 


(1,1) \ 


2 


i 

2 


1 





1 


i 

2 


+ 1 


1 


1 





1 

2 


1 


1 

2 


+ 1 


1 








1 

2 


1 


1 

2 








1 
















i 

2 





1 

2 








+1 



where Mi = ^(n 2 —n 3 ) , M 2 = 5(^3—^2)- The phases of the states \(A,0)nh} , |(0, A)i 
and |(cr, ^N^'Is) are those of Ref. [18|. The phase of the Clebsch-Gordan co- 
efficient is determined by forcing 



. // (A,0) (0,A) 
Slgn ,*I 2 



Icr, ol 



(2cr, o,0) \a 



(88) 



As always, it is convenient to rewrite Eq. ((87|) as 
\(a,a)N^I 3 ) = 5^|(i4,0)nJi)|(0,A)i/*/ 2 ) 



(A,0) (0,A) 



(o-,o) 



°/i Mi,/ 2 M 2 



(89) 



where ^ Ja M is the usual su(2) coupling coefficient and the reduced Clebsch 
Gordan/^ 2 ^ 

\ niii i/|7 2 



A = 2 and A = 1. 



^00^ ) does not depend on Ms. Tables arc provided for 
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Table 5: 2-atom case. Reduced 611(3) Clcsbsh-Gordan coefficients for (2, 0) <g> 
(0,2) -(0,0). 
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/(2,0) (0,2) 
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Table 6: 2-atom case. Reduced 6X1(3) Clcsbsh-Gordan coefficients for (2, 0) <£> 
(0,2) -,(1,1). 
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\n 1 I 1 v*I 2 
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Table 7: 2-atom case. Reduced ©11(3) Clesbsh-Gordan coefficients for (2, 0) <E> 
(0,2) -> (2,2). 
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C Final form of the density matrix for the non- 
degenerate case 

In this section we present the technical steps to obtain Eq. (f3l?]) from Eq.(|38p. 

First, multiply both sides of Eq.® by v[^ Q n mj,( t ) for fixed iV» and 
fixed J', integrate over the Sil(3)-invariant measure of Eq. (|7i|) and rearrange. 
This produces 



In this last expression, the sums over n and v are not independent but linked 
by Eq.dMl)- 

To complete the inversion, we use orthogonality of 611(3) CGs: 




(/w)0, NM J 



) 



(90) 




(A,0) (0,A) (fx,fx) 
nl v*V N(v)j> 



) 



( 



(A,0) (0,A) (fx,n) 
hi D*F NWJ' 



) 




(91) 



and rearrange the notation to finally yield Eq.(|39]). 
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